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lO ' Abstract 

The multimagnon continua of ID quantum spin systems possess several interesting singular 
\ features that may soon be accessible experimentally through inelastic neutron scattering. These 

' include cusps and composition discontinuities in the boundary envelopes of two-magnon contin- 

■ uum states and discontinuities in the density of states, "caustics" , on and within the continuum, 

which will appear as discontinuities in scattering intensity. In this note we discuss the general 
origins of these continuum features, and illustrate our results using the alternating Heisenberg 
antiferromagnetic chain and two-leg ladder as examples. 



1 Two-magnon states 



1.1 Introduction and definitions 

Recently the subject of two-magnon excitations of quasi-lD quantum spin systems has attracted 
considerable interest. Theorists have long predicted that some of these systems will possess bound 
states and there are now experimental indications of such bound states in the alternating chain 
material copper nitrate and in the spin ladder (Ca,La)i4Cu2404i P, ^. 

In addition to the bound modes there is a continuum of two-magnon states, which has been 
reported in inelastic neutron scattering from polycrystalline Sro.73Cu02 (assuming an alternating 
chain model) and copper nitrate single crystals 0. Although this continuum has been considered a 
rather uninteresting feature in comparison with the bound states, it actually possesses an interesting 
and rather complicated internal structure p, 0]. With high resolution inelastic neutron scattering 
experiments it should be possible to study these features of the two-magnon continuum in parallel 
with studies of the bound modes. To facilitate these experiments, here we discuss the nature of 
these interesting features of the continuum and show how they can be understood simply in terms 
of aspects of the one-magnon modes. 

The two-magnon continuum states in a gapped ID antiferromagnet are simply composed of two 
one-magnon excitations, since the interaction region is a zero matrix element for these unlocalized 
states. The energy of a two-magnon state composed of magnons with ki and k2 is therefore given by 

uj2{k = ki + = uj{ki) + uj{k2) (1) 

with the usual crystal momentum constraints on ki, ^2 and k, 

k = ki + k2 mod{2Ti) . (2) 

Thus the allowed states in the two-magnon continuum are given by the simple vector addition of 
Fig.l. 

For a fixed two-magnon k we can independently vary fci, keeping k2 = k — ki. This gives a range 
of allowed two-magnon energies for given k. Carrying out this "scan" in ki at each k gives the full 
set of states in the two-magnon continuum. Interesting physics questions regarding this continuum 
include the determination of the upper- and lower- limit boundary curves of this allowed region, and 
the composition of the two-magnon states along and within these boundaries. 

All this information is implicit in the one-magnon dispersion relation uj{k). In the following 
discussion unless otherwise specified we will assume that uj{k) is an even, periodic function of k 
with period 27r, monotonically increasing inside the range k = [0,7r], with a single inflection point 
at /cq. We also assume that uj{k) departs from its = minimum and k = ±7r maximum values 
quadratically in k, and that the inflection point lies within \k\ = [0,7r/2]). (This a;(/c) is abstracted 
from the simple alternating chain result, shown in Fig.l. Our discussion can easily be extended to 
more general cases.) 

1.2 Lower continuum boundary 

First we consider the lower boundary of the two-magnon continuum, which we call fig (^)- The A; = 
two-magnon mode clearly has a lower boundary at ci;2(0) = 2uj{0), since this satisfies the momentum 



constraint k = ki + k2 mod(27r) and is the global minimum of free two-magnon energies. For other 
k values we wish to solve 

= mink^^k2=k-ki{^{ki) +(^(^2)) ■ (3) 

Starting from the symmetric state {ki,k2) — {k/2,k/2), with energy uj2{k) — 2u{k/2), we can 
establish this lower boundary ^2 (^) by shifting these momenta to (/ci,fc2) = {k/2 — A,k/2 + A), 
and determining whether any nonzero, physically distinct A (0 < A < tt) gives a lower-energy two- 
magnon state. (See Fig. 2 for this construction.) The condition that the initial, symmetric state itself 
has the lowest energy is 

u{k/2 - A) + u;{k/2 + A) - 2u;{k/2) > V A , (4) 

which can be used as the definition of a "globally concave" function. (A function is locally concave 
if this is satisfied for infinitesimal A.) The lower boundary curve of the two-magnon continuum is 
given by 

n^{k) = 2cu{k/2) (5) 
if the one-magnon dispersion is globally concave. 

Both upper and lower boundaries will usually be extremal in the energy of two magnon states, 
so that they satisfy 

duj2{ki, A; - /ci) - - (u;'{ki) - u'{k - ki)) dki (6) 

or equivalently 

uj'{k/2-A)=uj'{k/2 + A). (7) 

The search for boundary curves of this extremal type can be considered a search for the minimum- 
or maximum-energy solutions A='=(A;) of this equation; once (7) is solved for A''=(A;), the boundary 
curves are then given by 

nt{k) = u{k/2 - A^{k)) +u{k/2 + A^{k)) . (8) 

Clearly one solution of (7) is A{k) = 0, giving Q2{k) = 2uj{k/2). This may however not give the global 
minimum or maximum two-magnon energy at k. We shall argue below that the lower boundary is 
in fact not given by A{k) = if u;{k/2) is convex. 

If the one-magnon dispersion relation uj{k) is not locally concave over the full physically indepen- 
dent range [— tt, tt], the concave region [— /cq, ^0] will be bounded by inflection points at k ~ ±/co. For 
a two-magnon state with momentum \k\ > 2kQ (which lies in the convex region of the one-magnon 
dispersion relation) the lower boundary of the two-magnon continuum is not given by 2uj{k/2), since 
an infinitesimal perturbation into an unsymmetric state (A;i,/c2) = {k/2 — A, k/2 + A) lowers the 
energy; 

lim^a^(A;/2- A)+u;(A;/2 + A) -2a;(A;/2) = AV'(A;/2)+0(A^) < . (9) 

The nature of the departure of f2^(A;) from the symmetric state as we enter the convex region 
depends on the behavior of u!{k) near the inflection point ko- If the departure from the symmetric 
state A{k) = is continuous, A(A;) will be small near k — 2ko, and we can expand (7) to flnd 



lim A(k) 

fc— >2feo_|_ 



-6m;(")(A;/2)]V2 ^^^^ 



w(^'')(A;/2) 



We can further simplify this by expanding in A; — 2/co and using cu'^ko) — 0, which gives 



hm A{k) = 

k—*2ko+ 



w(^^)(A;o) 



1/2 



(11) 



Thus if the fourth derivative c<j*^™)(A:o) is positive (c<j*^'"*)(/co) < at this inflection point), the minimum 
energy state (fci, ^2) = {k/2 — A{k),k/2 + A{k)) departs from the symmetric one {k/2,k/2) as a 
square root of momentum. 

The corresponding departure in energy of (^) from 2uj{k/2) is much less abrupt, and for small 
k — 2ko is 



lim Q^{k) -2uj{k/2) 

A;^2fco+ 



{k-2koy . 



(12) 



3 lw^'''\kof^ 

4 I. w^'^'^iko) 

This picccwisc continuous behavior and square-root departure in composition {ki, ^2) from sym- 
metric states on (A;) near 2ko is illustrated by the alternating chain in Fig. 3. In this model to 
leading order in a (see Appendix) one has /cq — 7r/2 = —a/2, w^^^^\ko) = —a/2 and w^^'"^ko) = 3q;^/4, 
so 



lim A(k) 

a— ►O 



TT — a 



(13) 



and 



Jm ^-^{k)-2u{k/2)^ _1 (A;-(vr-a))' . 



ct^O 



(14) 



The location of this composition discontinuity could be used as an independent experimental deter- 
mination of a that does not require following cj(A;) over a wide range of k. 

Between \k\ = 2ko and |A;| = tt the two-magnon lower boundary passes through states of increas- 



ingly unsymmetrical composition, reaching the limit ki = Q,\k2\ 
in Fig.3. 



TT at = TT. This is illustrated 



Were w^^^^ko) negative we would find a discontinuous jump in the value of ki on the boundary 
^2 (^) at A; = 2A;o. In either case we expect to see rapid variation in scattering intensities from states 
near this composition discontinuity. 



1.3 Upper continuum boundary 

The upper two-magnon continuum boundary, defined by 

^lt{k) = maxki^k2=k-ki (uj{ki) + uj{k2)^ , (15) 

can be found using a similar construction to Fig. 2. First note that the global maximum of energy 
of two-magnon states is 2a; (tt), ai k = 2n (equivalent to /c = or any k = 27im) and consists of 
a symmetric state of two magnons, ki = k2 = ti. As we move away from this highest-energy state 
at A; = 27r, we can form analogous symmetric states with total momentum 2t[ — k (equivalent to 
—k and degenerate with a k state) from two single-magnon states with A;i = A;2 = tt — A;/2, as 
shown in Fig.4. Evidently this state has energy 2uj{'k — A;/2). We can determine whether this is the 
maximum energy state at this k by considering a symmetrically displaced combination {ki,k2) — 
{n — k/2 — A, n — k/2 + A), analogous to Fig. 2. The energy difference is 



AE = u{n - k/2 + A) + uj{tt - k/2 - A) - 2uj{n - k/2) , 



(16) 



which is negative for all k in [0, it] (meaning that the symmetric state has the highest energy) if uj{k) 
is convex over the range [tt/2,7i]. If this is satisfied the upper boundary curve of the two-magnon 
continuum is given by 

nj{k) = 2uj{7v ~k/2) . (17) 

This is the case for our alternating chain example, as is shown in Fig. 4. It is also satisfied for the 
uniform S = 1/2 Heisenberg chain, for which Lj{k) = {tt/2)J \ sm{k)\\ since this is a convex function, 
(17) implies 

^]+(fc) = 2uj{ti - k/2) = 71 J I sin(A;/2)| , (18) 
which is the well-known upper boundary of the two-spinon continuum [Q. 

If instead we encounter an inflection point uj"{kQ) = in the interval [7r/2,7r], we will find a 
departure from 2u!{tt — k/2) along the upper boundary. As in Fig. 3 this will be accompanied by a 
change in composition from the symmetric state (fci, ^2) = (tt — A;/2, vr — /i;/2). 

Since we have assumed for our example that the inflection point k^ lies in the interval [0,7r/2], 
we do not encounter an inflection point, and the function VL2{k) = 2uj{7c — k/2) describes the entire 
upper boundary curve for k G [— 7r,7r]. This function is monotonically decreasing inside the range 
k = [0, vr], and has a negative slope at /c = vr, given our assumption that uj{k) has a positive slope at 
7r/2. Reflection symmetry of the dispersion relation about k = n then implies a cusp in the slope of 
the upper boundary curve cit h = IT, as shown in Fig. 4. 

In Fig. 5 we show the complete result for the boundary curves ^2 (^) of the two-magnon 

continuum for our alternating chain example. Note that the two-magnon bandwidth minimum at 
A; = vr is not zero; the value 

ly™" = fi^(7r) - Q^in) = 2uj{n/2) - w(0) - u{tt) , (19) 

is a^/4-|-a^/8-|-0(a;^) for the alternating chain. A zero bandwidth at /c = vr would follow for example 
from pure cos(A;) modulation in the one-magnon dispersion relation. 

1.4 More general boundary curves 

Motivated by the alternating chain as a prototypical gapped ID quantum spin system, we have 
assumed that the one-magnon dispersion relation is even, monotonically increasing inside [0, tt], and 
has a single inflection point at k^ inside [0,7r/2]. This led to an upper two-magnon continuum 
boundary with a cusp and a symmetric composition k\ = k2 everywhere, and a lower two-magnon 
boundary with a composition break at 2kQ. 

Although these results are valid in general for alternating chains, in other quantum spin systems 
we may encounter a single inflection point in [tt/2,tt], or a different number of inflection points. An 
example of more complicated behavior is provided by the two-leg ladder, which for sufficiently large 
a has a one-magnon dispersion relation uj{k) with db ixiiiiiixiuixi cit k = TT, a secondary minimum at 
k = 0, two inflection points in G [0, vr], and a maximum at an intermediate /c G [0, 7r/2]. 

We expect to find a symmetric state ki = k2 = k/2 at the lower boundary ^2 (/c) of the continuum 
when uj{k/2) is globally concave, and an unsymmetric two-magnon state on ^2 (/c) with an energy 
< 2uj{k/2) when uj{k/2) is a locally or globally convex function of k. With more than one inflec- 
tion point, boundary curves will typically be composed of pieces with symmetric and asymmetric 
composition. 



1.5 Density of states: caustics 



Our results for the composition of two-magnon states along the upper and lower continuum bound- 
aries imply the existence of discontinuities in the density of states on the continuum boundary and 
under certain conditions within the continuum. We refer to any such discontinuity in the density of 
states as a "caustic" . 

Consider the density of states in a system of N spins with unit cell b; the number of one-magnon 
states dn in an interval dk is 

dn = (^^^ dk = nodk . (20) 

and for two-magnon states in dkidk2 it is 

d^n = nl dkidk2 ■ (21) 

Since we observe two-magnon continuum states with specified total energy E = uj{ki) + uj{k2) + EQ 
and total momentum k = ki + k2, the number of states d^n in a region dkdE is more relevant to 
experiment. This density has a Jacobean factor, 

d'^n^ -dkdE = p(k,E) dkdE . (22) 
I u;'{ki) - uj'{k2) I 

Note that the condition for a divergent density of states p{k, E), which is a singular caustic, is 

cu'iki) - uj'{k2) , (23) 

which was exactly the condition (7) used to search for boundaries of the two-magnon continuum. 
Thus we expect to find a divergent density of states p{k, E) on the two-magnon boundary curves 
^2 (^) and nt{k). 

Unfortunately we cannot invert the {E, k) ^ (fci, ^2) relationship and solve (22) for p(/c, E) for a 
general one-magnon dispersion relation u}{k). We can however determine this density of two-magnon 
states in certain limits, including the case of states close to the boundaries. 

First we consider the symmetric lower boundary VL2{k) = 2uj{k/2). (Recall that 2uj{k/2) gives 
the lower boundary of the two-magnon continuum for a range of A;, as discussed in previous sections.) 
Perturbing the momenta away from this symmetric point into the continuum of asymmetric states, 
we can invert {E, k){ki, k2) infinitesimally and use (22) to find the density of states p{k,E) near 
2(jj{k/2). This gives 



n 



^J5V2)^^^'''^ = 2-|a;'W2)|V^ 



E - 2uj{k/2) 



-1/2 



(24) 



which evidently has a l/y/AE singularity as we approach 2u{k/2). This behavior is reminiscent of 
the van Hove singularity in the density of states as we approach a band edge, although that single 
particle density is instead proportional to l/\u!'{k)\. 

The asymmetric region of il2 (^) with ki < k2 also has a singular density of states, since u!'{ki) — 
uj'{k2) there as well. We may again solve for p{k, E) near the boundary using (22), with the result 



lim p{k, E) — 



nl 



E-.oj{ki)+oj{k2)'^^ ' ' I 2(a;"(A;i) +a;"(A;2)) |V2 



E — Lo{k\) — Lo{k2) 



-1/2 



(25) 



where {ki, k'i) — k/2 ^ A(A;). The symmetric boundary formula (24) is a special case of this more 
general result. 



A final special case is the density of states in E along the line k — 2ko; since u!"{ko) - 
singular case. Expanding (7) to O(A^) (the usual leading 0(A) term vanishes), we find 



lim p(k 



2ko, E) 



Tin 



1/4 



E - 2(jj{ko) 



-3/4 



this is a 



(26) 



.4|a;('^)(A;o) , 

Although we have specifically discussed the singular density of states found at the upper and 
lower boundaries ^lf(k), we note that there may be caustics specified by (7) within the continuum 
as well, and that these will also have square-root divergences in p{k, E) as we approach the singular 
line. 

We again use the alternating chain to illustrate this. In Fig.6 we show points in the two-magnon 
continuum, generated with a p{k, E) distribution. The divergences in p(/c, E) on the upper and lower 
boundaries, "border caustics", are evident (see also Fig. 10). Fig. 7 shows an enlargement including 
the asymmetric region k > 2/co; a caustic within the continuum due to the solution u!2{k) — 2uj{k/2) 
of (7) is apparent. 



1.6 Spin ladder 

As an independent example of a two-magnon continuum, in Fig.8 we show the one-magnon dispersion 
uj{k) and two-magnon continuum for a two-leg spin ladder, as in Fig.6 for the alternating chain. The 
coupling strengths are a = J\\/Ji_ = 0.3 and Jj^ = 1, and the analytic results used are summarized 
in the appendix. 

Although the ladder one-magnon ui{k) has the locations of maxima and minima translated by 
TT relative to the alternating chain, the two-magnon continua of the two systems are quahtatively 
remarkably similar for small a. This is because a shift in /c by 27r, for two magnons, is equivalent to 
zero. Close inspection reveals that the a-dependent effects are typically rather larger on the ladder 
than the alternating chain, presumably because the ladder dimers are coupled by twice as many links 
as the alternating-chain dimers. 

2 Higher multimagnon states 
2.1 Definitions 

Higher multimagnon continua show features similar to the cusps and discontinuities of two-magnon 
states, although these will presumably be of much less experimental interest due to the weaker 
coupling of higher-lying states to experimental probes. We shall nonetheless sketch the generalization 
of some of the simpler results for two-magnon systems to the n-magnon continua. 

The energy of an n-magnon state is given by 

UJn{k) = Uj{km) (27) 
m=l 

where 

n 

k^Y^^rn mod (27r) . (28) 

m=l 



Again assuming a gapped one-magnon dispersion relation with the general features discussed in 
the previous section (an even uj{k), monotonically increasing inside [0,7r], quadratic behavior in k 
near these endpoints, and a single inflection point ko), we can infer several features of the n-magnon 
continuum. 

2.2 N-magnon continuum boundaries 

Provided that the individual magnon momenta k/n are below the inflection point ko, {k < nko V A; e 
[0, tt] ) , the lower boundary fl~ {k) of the n-magnon continuum is given by the energy of the symmetric 
state (/ci, k2. ■ ■ kn) — {k/n, k/n .. .k/n), 

n-{k) = nw{k/n) (29) 

This implies that for sufficiently large n {n > n/ko) all the lower-boundary curves {Q~{k)} possess 
cusps at A; = TT. Since the departure of u{k) from a;(0) for small k is assumed to be of the form 
u!{k) = a;(0) -|- C2k^ + 0(A;^), the lower boundary curve will approach 

\mi^n-(k) = nw(0) + (c2/n)k'^ + 0(n-^) . (30) 

Thus ^n{^) asymptotically approaches the line nu!{0), with a residual 0{n~^) parabolic component 
oc k'^/n. 

The upper boundary f^^(/c) differs for even and odd magnon number. For an even number of 
magnons the upper boundary maximum energy nuj{7r) occurs at /c = nvr = 0. Again assuming 
quadratic behavior near the extremum, u{k) = a;(7r) — 02(71 — kY + 0((7r — A;)^), the symmetric 
solution {ki, /c2 • • • kn) = (tt — k/n, n — k/n . . .tt — k/n) with energy 

n+{k) ^nu;{7r-k/n) (31) 

gives the upper boundary curve provided that we do not encounter an inflection point, n — k/n > kg. 
Once again this will be satisfied for sufficiently large n, and all these even-n upper boundary curves 
have a cusp at A; = tt. The behavior of this upper boundary for large n approaches a constant with 
a negative parabolic component of 0{n~^), 

n+{k) = nw{7r) - {c2/n)k^ + 0(n-=^) , (32) 

analogous to the result for the lower boundary. The upper boundary for odd magnon number n 
differs in that the maximum energy state has total momentum k — tt rather than A; = 0, and a cusp 
at A; = rather than k — tt. 

2.3 Caustics 

The higher multimagnon continua also possess discontinuities in their densities of states {p{k,E)}, 
although these are qualitatively different from the results we have shown for the two-magnon case. 
This is largely due to the fact that the {ki,k2) —>■ {k,E) mapping gives a simple Jacobean (22) 
which leads to caustics where it is singular, whereas the n > 3 multimagnon continua project a fiat 
distribution in {ki, k2, ■ ■ ■ kn) onto a smaller number of variables (A;, E). This involves an integration 
over n — 2 variables, which smoothes the singularities seen in the two-magnon case. 



As examples of these higher continua, in Fig. 9 we show the three-magnon continuum for our 
alternating chain example, and Fig. 10 shows the density of states encountered in a fixed- A; slice 
(at k = 7r/2) through the two-, three-, and four-magnon continua. (This figure is a histogram of 
points falling within 0.4957r < k < O.SOStt, with an energy binning of 6E = 0.01, normalized as 
f{k,E) — N{k,E)/SE 5k Ntot- The binned points were selected from an initial flat distribution, 
— TT < ki,k2 . . .kn < oi Ntot. — points.) 

The two-magnon continuum in Figs. 9, 10 clearly shows the l/y/AE singular border caustics in 
the density of states, as given by (24). In contrast, the three-magnon continuum density p{k, E) has 
border caustics with only finite discontinuities in the density of states, and the four-magnon p{k, E) 
is continuous on the boundaries. 

3 Summary and Conclusions 

We have considered the multimagnon continua of ID quantum spin systems, and have illustrated 
our results using the Heisenberg alternating chain and two-leg spin ladder as examples. 

The boundaries of the two-magnon continua that result from one-magnon dispersion relations 
similar to the alternating chain case are considered in detail, and general results are derived for 
these and higher multimagnon boundary curves. It is shown that these curves exhibit cusps and 

discontinuous changes in composition under certain conditions. The density of states p{k, E) in the 
continuum is also considered, and it is noted that p(k, E) can possess discontinuities and divergences 
on "caustic" lines, both on the continuum boundaries and within the continuum. 

We anticipate that it should be possible to identify these features of two-magnon continua through 
high-resolution inelastic neutron scattering experiments on candidate alternating chain and ladder 
materials, and that precise studies of these features may allow determinations of the parameters of 
spin Hamiltonians through the study of relatively small ranges of {k, E) space. 

4 Acknowledgements 

It is a pleasure to thank B. Lake, G.l. Meijer, G. Miiller, S.E. Nagler, T. Papcnbrock, J. Ricra, R.R.P. 
Singh, D.A. Tennant and G.S. Uhrig for useful discussions and communications. This research was 
sponsored by the Laboratory Directed Research and Development Program of Oak Ridge National 
Laboratory (ORNL), managed by UT-Battelle, LLC for the U.S. Department of Energy under Con- 
tract No. DE-AC05-00OR22725, and by the Department of Physics and Astronomy, the Neutron 
Sciences Consortium, and the Chemical Physics Program of the University of Tennessee. 



Appendix: results for alternating chains and ladders 



In this appendix we collect some results for the alternating Heisenberg chain and two-leg ladder, 
which we used to illustrate aspects of multimagnon continuum states in the paper. 



The alternating chain Hamiltonian is defined by 

N/2 

H = S2i-i ■ S'2i + a S2i ■ S2i+i , (Al) 

i=l 

where we assume cyclic boundary conditions. The corresponding one-magnon dispersion relation 
uj{k) can be written as a Fourier series, with a-dependent coefficients; 

oo 

uj{k) =Y,ai{a) cos{ik) , (A2) 

where we have implicitly set the strength J and unit cell b to unity. The Fourier coefficients may 
be evaluated as a power series in a using a strong-couphng expansion. The O(a^) series of Ref.Q, 
which we use to generate the numerical results given in this paper, are 
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(A3) 
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We may use these coefficients to develop series for several of the quantities discussed here. In 
particular, the inflection point of uj{k) is given by 

ko = ---a~-a^ -}}^a^ + 0{a'>) (A4) 
" 2 2 96 576 ^ ' ^ ' 

and the energy (^2(2^0) at the "break point" 2kQ (see Fig. 3) on the lower boundary of the two-magnon 
continuum is given by 

1 3 11 31 1 

Q- (2ko) = 2uj(ko) = 2--a^-—a^ - - — + O(a^) . (A5) 

2 V 0; V 0; 2 32 64 1024 ^ ^ ^ ^ ' 

The derivatives of the one-magnon energy ui{k) at k^ determine the nature of the asymmetric two- 



magnon lower boundary state for k > 2ko. To 0{a^) these are 
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We also give the corresponding results for the two-leg spin ladder (to one higher order in a) . The 
A'"-rung ladder Hamiltonian is given by 

H — J±'^ Si^i ■ Si^2 + OL {Si^i • Si+i^i + Si^2 • 'S'i+1,2) , (A7) 

i=l 

and the u;{k) one-magnon Fourier coefficients to 0{a^) are: 
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The inflection point ko is: 

A;o = I + « + + ^Q;'-^Qi^-^Qi' + 0(Q;^) (A9) 

and the lower boundary energy ^2(2^0) the break point 2kQ is: 

n^{2ko) = 2u;{ko) ^2-^a' - ^a' - ^a' + ^a' + 0{a') . (AlO) 



Finally, the derivatives of a; (A;) at ko are: 



u'iko) = -a +la^ -^a^ -Ti«^ -fi"^ 

cu"(A;o) = 

a;(-)()to) = +a -la' +§a' a' 

u;(-)(A;o) = +3^2 ^3^3 +|^4 ^m^5 +^^6 



(All) 
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Figure 1. The alternating chain one-magnon dispersion relation uj{k) for J = 1 and a = 0.3, 

showing a symmetric two-magnon state. 




Figure 2. Illustrating the construction of an asymmetric two-magnon state (4). 




Figure 3. The k = 2ko transition from symmetric (fci = ^2) to asymmetric {ki < /C2) states on the 

lower boundary (/c) of the two-magnon continuum. (Ahernating chain, J = 1 and a = 0.3; 
2fco/7r = 0.898, fi^(2A;o) = 1.950.) The composition break in + /C2) is shown at lower right. 
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Figure 4. Construction of the two-magnon upper boundary i72^(A;). (Alternating chain, J = 1 and 

a = 0.3.) 
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Figure 5. The complete two-magnon continuum boundary fi2^(fc) and ^2 (A;), showing the 
composition break at = 2fco- (Alternating chain, J = 1 and a = 0.3.) 




Figure 6. The density of two-magnon states within the continuum. (Alternating chain, J = 1 and 

a = 0.3; flat distribution in ki and /c2 of 2^^ states.) 
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Figure 7. An enlargement of the alternating chain continuum of Fig.6, showing a discontinuity in 
the density of states p{k,E) (a "caustic") within the two-magnon continuum. The caustic is the 
line 2uj{k/2) and its reflection about k = n. The composition break at A; = 2fco where this curve 
departs from the lower continuum boundary ^2 (A;) is also indicated. (Flat distribution in ki and k2 

of 2^^ states in the region displayed.) 
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Figure 8. The one-magnon uj{k) and density of states within the two-magnon continuum for a 
two-leg spin ladder. {J± = 1 and a = 0.3; fiat distribution in ki and k2 of 2^^ states.) 
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Figure 9. Two- and three-magnon continua of the alternating chain. ( J = 1 and a = 0.3; flat 
distribution in ki and k2 of 2^^ states in each case.) 




Figure 10. The density of states p{k,E) encountered in a fixed-k slice {k = 7r/2) through the two-, 
three- and four-magnon continua. (Alternating chain, J = 1 and a = 0.3.) 



